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VOLUMETRIC LAMINATED SAND ANALYSIS 

[0001] This application claims the benefit of U.S. Provisional Application Number 
60/265,659 filed January 31, 2001. 

BACKGROUND OF THE INVENTION 

[0002] This invention relates generally to methods for analyzing data obtained from 
well logs to determine the hydrocarbon pore volume of a subsurface geological 
formation. More specifically, this invention relates to an analysis of well log data 
taken in a reservoir having thinly interbedded sandstone and shale layers. 

[0003] An estimate of the hydrocarbon pore volume (HPV) of a subsurface 
geological formation has been obtained by the analysis of well log data obtained from 
an exploration well penetrating the formation. A well log generally involves the 
recording of a series of measurements obtained while an instrument traverses the 
formation through a wellbore. Such measurements are recorded in correlation with 
uniform depth increments. The well logging instruments produce measurements that 
represent an average of the measured rock property over a finite length of the wellbore 
which typically is about three to six feet. The various types of well logging 
instruments commonly used make measurements based on electrical, acoustic, and 
nuclear properties of the formation. 

[0004] In reservoirs where sand and shale are interbedded at bed thicknesses of two 
feet or less, the vertical resolution limits of the well logging instruments result in HPV 
estimates being too low. Accordingly, there is a significant risk that economically 
significant hydrocarbon reserves may be overlooked. Also, if an induction-logging 
instrument is used to make electrical resistivity measurements in a thinly bedded 
reservoir, the resistivity measurement becomes less sensitive to the hydrocarbon 
saturation in the sandstone layers and more sensitive to the proportion of shale. As a 
consequence, there is less certainty in the estimate of HPV in thinly bedded reservoirs. 
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[0005] Heretofore, two general approaches have been developed to seek more 
accurate estimates of HPV from well log data obtained in thinly bedded reservoirs. 
One approach is referred to as "high-resolution," and the other approach is referred to 
as "low-resolution." The high-resolution approach seeks to enhance the vertical 
resolution of the logging instruments used. The vertical resolution of a logging 
instrument refers to the thinnest bed in which a true reading can be obtained. The 
low-resolution approach seeks to determine the relationship between HPV and 
averaged well log data over a gross formation interval (e.g., 10 feet or more in 
thickness) so as to obviate having to resolve accurate log data values in individual thin 
beds. 

[0006] The essential element of the low-resolution approach is in establishing the 
relationship within overall interval conductivity of sand layer conductivity and shale 
layer conductivity. In the simplest case, this relationship can be expressed as: 
C int = h sand x C sand + x C shale , where C int , C sand , and C shale are the interval, sand, and 
shale conductivities, respectively; and h sand and h shale are, respectively, the fractional 
sand and shale thickness within the interval of interest, with h sand + h shale = 1. The 
interval may be the resolution interval of the resistivity log or may be a larger 
reservoir interval. Similar equations, well known to those skilled in the art, exist for 
associated logs such as density, neutron porosity, and gamma ray, which may be used 
in the estimation of interval hydrocarbon pore volume. Illustrative of the prior art 
implementations of the low resolution approach are Hagiwara, U.S. Patent No. 
4,739,255, and Frans G. Van den Berg et al., "Sandwich: Log Evaluation In Laminate 
Shaly Sands," SPWLA 37th Annual Logging Symposium, June 16-19, 1996. 
[0007] Also, even though thinly bedded sandstone formations have historically 
been called shaly sands in the petrophysical literature, standard shaly sand log 
analysis techniques do not correctly account for the effect of thin bedding on log 
responses. The standard shaly sand log analysis techniques were developed to address 
the effects of dispersed clay in sandstone rather than macro scopically interbedded 
sands and shales. The electrical effects of these two modes of clay or shale 
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distribution are significantly different, and therefore, different log analysis techniques 
are required to correctly analyze each mode. 

SUMMARY OF THE INVENTION 

[0008] The present invention provides a low-resolution approach for thin bed 
analysis, which will be referred to herein as Volumetric Laminated Sand Analysis 
("VLSA"). The approach of the present invention provides for an analysis of standard 
well log data obtained in thinly bedded reservoirs to obtain improved estimates of 
hydrocarbon pore volume. The analysis is applied on the basis that the reservoir 
formation is constituted by a sequence of approximately parallel, planar beds that are 
classified into bed types and wherein each bed type can be characterized as to 
porosity, capillary pressure behavior (i.e., water saturation), and aspect ratio. In 
applying the analysis, an initial geological formation model (including bed types, bed 
properties, and bed frequencies) is produced for a subject well and consistent with its 
measured log data. Thereafter, a Monte Carlo inversion is performed to find the 
resultant probability distribution of HPV. 

[0009] In accordance with the VLSA methodology, the initial earth model is 
obtained by specifying the analysis interval of the thinly bedded reservoir and 
specifying a set of bed types (e.g., sandstone beds and shale beds) believed to 
constitute the analysis interval. For each bed type, estimates are made of the average 
values of the relevant petrophysical properties (e.g., porosity, water saturation, and 
thickness within the interval), aspect ratio, and any other parameter as may be desired 
based on core data or regional knowledge. Estimates are made of the frequency of 
occurrence of the different bed types, which may be estimated from conventional or 
high-resolution logs or from core image analysis. These estimates are combined, 
using appropriate theoretical models, to produce estimates for the average log values 
over the analysis interval. Averages of the actual log measurements are computed 
over the analysis interval. The estimated log measurements are compared to the 
actual log measurements to ascertain consistency. The model continues to be refined 
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by modifying the parameters until a model is found which is consistent with the 
measured log data. 

[0010] In the Monte Carlo inversion, a random earth model is generated consisting 
of independent random samples from the estimated probability distribution for all of 
the petrophysical parameters describing the initial model. Estimated log responses are 
calculated for this random earth model. The estimated log responses are compared to 
the measured log responses to determine consistency. If estimated and measured logs 
are consistent, the random earth model is retained. The Monte Carlo model continues 
to be refined by repeating the procedure over a predetermined number of trial cases. 
From the various retained cases, distribution statistics are computed for interval HPV 
and related parameters (e.g., sand porosity, sand water saturation, net-to-gross ratio). 
The distribution statistics represent the expected value and uncertainty in HPV 
associated with the average measured log values over the thinly bedded reservoir 
interval. 

[0011] Among the advantages of the VLSA methodology over the high-resolution 
approach and prior low-resolution approaches is that thin-bed effects on all required 
logs, and not just resistivity logs, can be taken into account. Also, the VLSA 
methodology accounts for the effect on the resistivity log of bedding units or 
inclusions with large aspect ratios (e.g., discontinuous beds or shale clasts). Further, 
in contrast to the high-resolution approach wherein new random uncertainty is 
introduced at each log depth, the VLSA methodology reduces such random 
uncertainty by consolidating data from multiple depths. 

[0012] Among the additional advantages of the VLSA methodology over prior low- 
resolution approaches is that deterministic, potentially unstable methods for 
calculating sand conductivity are not used. An example of such a potentially unstable 
method is the algebraic solution of the equation for interval conductivity set forth 
above. Instead, a completely general solution technique is used that provides stable 
solutions and uncertainty assessments over the complete range of shale fractions. The 
general solution technique extracts information on shale fraction from all available 
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logs, while also integrating any a priori information on shale fraction such as that 
obtained from a high-resolution image log. This avoids use of a pre-defined and 
limited set of log measurements (e.g., gamma ray, density, neutron porosity) to 
estimate the critical shale fraction parameter. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0013] Figure 1 illustrates sand-shale formations having two end-member bed types 
(sandstone and shale) and bed aspect ratios varying from 0.0 to 0.5; 

[0014] Figure 2 illustrates a borehole penetrating a series of approximately parallel, 
planar beds; 

[0015] Figure 3 illustrates the principal steps of the VLSA methodology; 

[0016] Figure 4 illustrates the Monte Carlo inversion procedure used in the VLSA 
methodology; 

[0017] Figure 5 is a plot of synthetic gamma ray, induction log conductivity, 
neutron porosity, and bulk density logs for a synthetic formation used to illustrate the 
VLSA methodology; 

[0018] Figure 6 is a tabulation of the formation data that summarizes the known 
petrophysical properties of the synthetic formation used to illustrate the VLSA 
methodology; 

[0019] Figure 7 is the model setup worksheet for an EXCEL® workbook 
implementation of the VLSA methodology, as applied to the illustrative example that 
is the subject of Figures 5 and 6; 

[0020] Figure 8 illustrates the model setup worksheet, applicable to the illustrative 
example of Figures 5, 6, and 7, for specifying the a priori distributions of the VLSA 
model input parameters in the form of a bounded normal distribution; 

[0021] Figure 9 illustrates the form of the distribution of an interval-average 
parameter between its minimum and its maximum in relation to total water saturation 
of high-quality sand; 
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[0022] Figure 10 illustrates a worksheet for viewing the parameters and calculated 
values from any one successful Monte Carlo trial; and 

[0023] Figure 1 1 illustrates the statistics generated by the Monte Carlo inversion 
analysis on the distributions of feasible solutions for hydrocarbon pore volume, sand 
fraction, total porosity, and total water saturation, and wherein the shaded histograms 
represent the distribution of feasible solutions while the empty histograms represent 
the a priori distributions for each quantity. 

DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT 

Introduction 

[0024] Prior to describing the VLSA methodology, certain basic principles of log 
analysis will be described. Petrophysical properties are intrinsic properties of a 
volume of rock, which have well-defined values independent of measurement 
techniques. These properties may be isotropic (i.e., independent of orientation) or 
anisotropic (value depends on orientation within the rock volume). Examples of 
isotropic petrophysical properties are total porosity and water saturation. Examples of 
anisotropic properties are conductivity, acoustic velocity, and permeability. In 
idealized cases, the isotropic petrophysical properties of a thin-bedded rock volume 
can be inferred from the corresponding properties of the individual beds via simple 
averaging relationships known to those skilled in the art. For anisotropic properties, 
such relationships are more complex. 

[0025] Log responses are estimates of petrophysical properties as measured by well 
logging tools. These responses depend in a complex way on the geometries of the 
tool and the borehole environment. In most cases, however, thin-bed log response can 
be adequately approximated using the idealized physical averaging relationship that 
applies to the underlying petrophysical property. All of the thin-bed log analysis 
techniques based on conventional logs use these simplified averaging relationships. 
The types of petrophysical averaging relationships include simple volume-weighted 
averages of isotropic properties (e.g., total porosity, total water saturation, bulk 



Express Mail No.: ET608077262US 



Date Filed: 01/09/02 



ATTORNEY DOCKET NO.: PM 98.086 



density, and matrix density); mass-weighted averages of isotropic properties (e.g., 
gamma ray); and orientation dependent averages of anisotropic properties (e.g., 
conductivity). 

[0026] Certain terminology also applies to the expression of relationships between 
petrophysical properties and the composite properties of thinly bedded formations 
measured at a scale greater than or equal to log resolution. These definitions refer to a 
fixed volume "V" of the thin-bed formation. The fixed volume V may represent an 
interval of core, a cylinder around the wellbore of given radius and height, or the 
investigation volume of a given logging tool. A thin-bed property "P" is the value of 
a petrophysical property or log response measured as an average over the total fixed 
volume. An end member property "P" is the value of a petrophysical property 
measured over that fraction of the fixed volume composed of a given end member 
rock type (the i-th rock type). The volume fraction "V" of the i-th end member is the 
ratio of the volume of that end member to the fixed total volume (pore space 
included). If there are N end members in a given sample, the volume fractions satisfy 
the equation: Vj + V 2 + . . .+ V N = 1 . The mass fraction "W" of the i-th end member is 
the ratio of the mass of that end member to the total mass of the fixed volume (pore 
space included). If there are N end members in a given sample, the mass fractions 
satisfy the equation: Wj + W 2 + . . .+ W N = 1 . 

[0027] Certain equations relating to thin-bed log responses of interest include: 
Bulk Density RHOB = VjRHOB^ V 2 RHOB 2 

Gamma Ray GR = W,GR, + W 2 GR 2 

Conductivity (parallel) C P = VjC, + V 2 C 2 
Conductivity (series) 1/C S = Vj/C, + V 2 /C 2 
[0028] Fig. 1 shows thinly bedded formations with two end member bed types 
(sand 10 and shale 12). Fig. 1 illustrates beds with a range of aspect ratios, as 
indicated. Here the aspect ratio is defined as the ratio of bed thickness to bed width, 
and the beds are idealized as oblate ellipsoidal bodies. An aspect ratio of 0.0 
represents parallel, planar beds and an aspect ratio of 1.0 represents approximately 
spherical inclusions. The average vertical and horizontal conductivity of formations 
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approximated by this idealized structure, which have positive (i.e., non-zero) aspect 
ratios, can be estimated using effective medium theory (EMT), as expressed in the 
following set of implicit equations: 

k , where 
f _ v k /R, 



Z( V i/ R i,v) 
i 

R k ,v=i+kv(^K-i) 5 

kv=(l + e k a Xe k -tan- 1 (e ll )Vei! ^ 



k , where 

„ _ "k/R>,h 

i 

R k , h =1 + L kh (a k /a p -l) 5and 
L, h =(l-L k>v )/2 

Inputs to these equations are the following terms: 

ot k = l/(aspect ratio of bed type k), 

a k = conductivity of bed type k, and 

v k = volume fraction of bed type k. 
Outputs from these equations, which can be obtained by a recursive computation, are 
the conductivities parallel and transverse to the major axes of the ellipsoids, cp and 
as, respectively. In the case of parallel planar beds (i.e., beds having an aspect ratio 
of zero), these correspond to the ordinary parallel and series conductivities and the 
EMT equations reduce to the equations given above for parallel and series 
conductivity of parallel planar beds. 
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[0029] Figure 2 illustrates parallel planar sand and shale beds 10 and 12 dipping 
with respect to a borehole 14, with relative dip angle theta (6). In this situation, for 
beds of any aspect ratio, the average conductivity measured by an induction log (mid) 
is: 

a ild =[o-^cos 2 (9)+a p a s sin 2 (e)] 1/2 
Volumetric Laminated Sand Analysis (VLSA) 

[0030] The VLSA methodology provides for analyzing standard well log data from 
thinly bedded reservoirs to obtain unbiased estimates of hydrocarbon pore volume. 
The method provides for incorporating all known information about the reservoir's 
petrophysical properties (e.g. porosity, water saturation, net sand thickness) into the 
log analysis. The method also provides for incorporating information on the average 
aspect ratios of the sand and shale beds. In addition, the method estimates the 
inherent uncertainty in the estimated hydrocarbon pore volume and in the associated 
reservoir parameters (net-to-gross ratio, porosity, and water saturation). 

[0031] Figure 3 illustrates the principal steps of the VLSA methodology in its most 
generic form. The steps outlined in Figure 3 are described in more detail below. 

[0032] Step 1: Formulate a priori petrophysical model of thin-bed formation. 

a. In a subject well, identify the top and bottom of the thinly bedded 
reservoir interval that is to be analyzed. 

b. Identify a set of bed-types (an earth model) which are assumed to 
comprise the thinly bedded reservoir interval. 

c. For each bed-type, establish a priori probability distributions 
describing the possible ranges of porosity, water saturation, relative thickness 
of the bed-type within the composite interval, aspect ratio of the bed type, and 
such other parameters as may be required to calculate the theoretical response 
of each log measurement to that bed-type. In establishing these a priori 
distributions, all available auxiliary data may be employed (including, but not 
limited to, core measurements of porosity, water saturation, and capillary 
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pressure, and estimates of bed-type relative thickness based on high-resolution 
image logs, core descriptions, digital analysis of core photographs, or results 
of standard log analysis in adjacent thickly bedded intervals). 

d. The petrophysical model space, M, is the vector space consisting of all 
possible combinations of values of the parameters identified in step c. If the 
number of parameters is N, then M is an N-dimensional space. 

[0033] Step 2: Obtain log data at chosen resolution. 

a. Select the well log measurements that are to be used to constrain the 
determination of hydrocarbon pore volume for the subject interval. This set of 
measurements must include a deep-reading conductivity log and may include 
any or all of the other log measurements normally run for this purpose. The 
log-data space, D, is the vector space of all possible combinations of these log 
values. If there are P logs, then D is a P-dimensional space 

b. Select one of two resolution modes (i or ii below) for the analysis and 
apply linear filters to all logs to bring them into consistency with the selected 
resolution mode. 

i) Total-interval resolution mode. In this mode, a single solution 
is obtained for the hydrocarbon pore volume of the identified reservoir 
interval. Each log is "filtered" by calculating its average value over the 
interval. The set of average log values is represented by one vector d 
in the log-data space, D. 

ii) Resolution adjusted to match the lowest-resolution log. In this 
case, appropriate linear filters are applied to each log so the filtered 
logs are all resolution-matched to the lowest-resolution log. Each 
logged depth is represented by a separate vector in the log-data space, 
D. 
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[0034] Step 3: Formulate log-response forward model. For each selected well 
log type, formulate a theoretical predictor ("forward model") that predicts the 
measured log value for each point in the petrophysical model space, M. In 
general, this formulation proceeds in two steps: a first equation or algorithm 
predicts the log value within each of the bed types; and a second equation or 
algorithm predicts the log value resulting from the interlayering of all the bed 
types defined in the petrophysical model. As an example, the conductivity log 
is discussed below. Other log types would proceed in a similar manner. 

a. Conductivity within each bed type . There are several well-known 
equations that may be used to predict conductivity of a homogeneous, 
electrically isotropic rock with known petrophysical parameters. These 
include the Archie, Waxman-Smits, or Dual Water equations, all of which are 
familiar to those skilled in the art of log analysis. 

b. Conductivity of interlavered thin beds . Equations describing the 
average induction-log conductivity as a function of the conductivities and 
relative frequencies of the bedding units are described above. These equations 
account for the effect of beds dipping relative to the wellbore, and the effect of 
beds with finite aspect ratio. 

[0035] Step 4: Define measure of fit between measured and modeled logs. 
Define a function which measures the closeness of fit between two vectors (a 
measured point dj and a predicted point d 2 ) in the log-data space, D. The 
measure function can be written e(dj,d 2 ). An example of such a function is the 
joint probability index: 

e(d 1 ,d 2 ) = nexp(-0.5((d 1ik - d 2 , k ) /std k ) 2 ) 

k 

In this equation, d 1Jc and d 2Jc are the k-th of P well log values belonging to the 
measured and predicted log vectors d t and d 2 respectively, and std k is the 
corresponding standard deviation which represents the combined measurement 
and model error for the k-th log. The joint probability index is proportional to 
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the probability that d x equals d 2 , assuming that the combined measurement and 
model errors are distributed normally for each log and are independent 
between logs. 

[0036J Step 5: Perform Monte Carlo inversion of forward model. Inversion is 
performed for each measured data-point d in the log-data space (see step 2). 
Inversion consists of the following steps, which are illustrated in Figure 4. 

a. Using the a priori distributions 16 established in step lc, the Monte 
Carlo trial generator 18 generates a random vector, m, in model space, M. 

b. Apply the log forward model equations 20 to the model vector m, 
obtaining a predicted vector d(m) in log data space, D. 

c. Compute the measure of fit e(d, d(m)) between measured and 
predicted log vectors (reference numeral 22). 

d. Keep (reference numeral 24) or reject (reference numeral 26) the model 
vector m using a rule based on the measure of fit e(d, d(m)). As an example 
of such a rule, if e(d, d(m)) is defined as in the equation in step 4 above, then 
m may be recorded in the solution dataset, S, with a frequency N ranging from 
0 to 10 given by: 

N = round(10-e(d, d(m))) 

e. Steps a-d are repeated for a prescribed number of trials, or until the 
solution dataset S reaches a prescribed minimum size. 

[0037] Step 6: Analyze solution dataset obtained from inversion. The solution 
dataset, S, represents a random sample drawn from the a posteriori probability 
distribution of petrophysical models; that is, the distribution of models that is 
possible after taking into account the constraints imposed by the well log data. 
Thus S can be used to calculate standard statistical measures of central 
tendency and dispersion of the hydrocarbon pore volume and its components 
(net-to-gross, porosity, and water saturation). These measures of dispersion 
represent the uncertainty inherent in estimating hydrocarbon pore volume from 
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well log data in thinly bedded reservoirs. As an example, the median (50th 
percentile statistic) may be used as a measure of the most likely hydrocarbon 
pore volume, while the 10th percentile and 90th percentile may be used to 
describe the range of possible values. 

[0038] The solution dataset, S, may also be used to examine correlation among 
model variables to evaluate the sensitivity of hydrocarbon pore volume (and related 
parameters) to other model parameters. As an example, it is often found that the 
hydrocarbon pore volume depends sensitively on the assumed value of the shale 
conductivity. Examination of this sensitivity in the sample dataset can be a useful 
a posteriori check on the validity of assumptions made in developing the a priori 
probability distributions (step lc). 

Synthetic Data Interval-Average VLSA Example 

[0039] A. Model Set-up. In Figure 5, synthetic gamma ray (GR) 28, induction 
log conductivity (C(induction)) 30, bulk density (RHOB) 32, and neutron porosity 
(NPHI) 34 logs for a synthetic formation are shown. The synthetic formation 
comprises three, interlayered bed types: a high-quality sandstone 36, a low-quality 
sandstone 38, and shale 40. The shaded areas in the left-hand column indicate the 
distribution of high-quality and low-quality sandstone and shale at 0.5 foot resolution. 
The interval to be analyzed is interval 42, as shown in the second column. The data 
set tabulated in Figure 6 shows the petrophysical properties of these bed types, 
including average total porosity water saturation (SJ, cementation exponent (m), 
saturation exponent (n), aspect ratio, volume fraction, brine conductivity (C w ), angle 
(0), interval thickness (feet), and hydrocarbon pore volume (feet). 

[0040] In Figure 7, a model of the formation based on the illustrative example of 
Figures 5 and 6 is presented in a spreadsheet implementation of the average-interval 
approach of the VLSA methodology. The model is a starting point for the Monte 
Carlo analysis. The preferred spreadsheet implementation uses EXCEL® version 5.0 
or higher and has the EXCEL SOLVER® installed. The worksheet data includes input 
petrophysical parameters, input log means and standard deviations, calculated 
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petrophysical averages, and calculated log relative errors and cumulative errors. The 
bed type fractions and calculated averages represent a base case solution to the thin- 
bed problem, matching all three log inputs with relative errors less than 1.0 in 
absolute value. In the model set up, the parameters are as follows: 

PHIT = porosity 

= total water saturation 
S wb - bound water saturation 
m = Archie cementation exponent 
n = Archie saturation exponent 
rhoma = matrix density 
rhof = fluid density 
GR = average gamma ray log values 
RHOB = bulk density (p b ) 
C t = conductivity 

Aspect = aspect ratio of the bedding units 

frac = bed type fraction 

HPV = hydrocarbon pore volume 

= free- water conductivity 
C wb = bound-water conductivity 

angle = the relative angle between bedding planes and borehole 

h = thickness of the analysis interval 
[0041] The log responses for each bed type are computed using standard 
petrophysical relationships as follows: 

Bulk density RHOB = (PHIT x rhof) + [(1-PHIT) rhoma] 

Conductivity C t = [C^ + S wb (C wb - C^/ SJ x PHIT m SJ 1 
The thin-bed model averages are obtained as described above. Hydrocarbon Pore 
Volume is computed as: HPV (feet) = h(feet) x PHIT x (1- SJ. HPV is the main 
parameter that is being sought through the VLSA methodology. Log relative errors 
are computed as (Model average - Input log mean)/Std. Dev. The mean squared error 
(MSE) is the mean of the squared log relative errors. The maximum absolute error 
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(MAXABS) is the maximum of the absolute values of the log relative errors. If 
MAXABS is less than one, then all the computed model averages match the input log 
means to within one standard deviation. 

[0042] The objective to finding a base-case solution is to adjust bed-type input 
parameters and volume fractions to achieve a model where MAXABS is less than one. 
That is, all the predicted logs fit the input logs to within the specified standard 
deviation tolerances. Bed-type parameters and volume fractions can be adjusted 
manually to reduce MAXABS and seek a base-case solution. Such adjustments 
should consider the relative confidence in the different parameters based upon the 
source. For example, if the average porosity for high-quality sand (HiQ Sand) was 
obtained from a large core plug database that showed relatively little variation in 
porosity, then little or no adjustment should be made to that parameter. The focus 
should be on adjusting those parameters that are known with the least certainty. 
During the iterative process, EXCEL SOLVER® can be invoked to vary the volume 
fractions to minimize MSE. If, however, the volume fractions are known with a high 
degree of confidence, such as from core image analysis, then they should be entered 
manually and not adjusted. Both manual adjustment and use of EXCEL SOLVER® 
can be undertaken until an acceptable solution is found. The base-case solution forms 
the basis for the subsequent Monte Carlo analysis. 

[0043] The base-case solution for the illustrative example is shown in Figure 7. 
This solution was found by using the EXCEL SOLVER® to adjust the volume 
fractions. It can be seen that this base-case solution comes very close to duplicating 
the correct formation parameters, which are illustrated in Figure 6. 
[0044] B. Monte Carlo Analysis. Monte Carlo analysis is a statistical method 
that yields probability distributions for unknown parameters. Further details 
concerning Monte Carlo methods can be obtained by reference to Tarantola, A., 
"Inverse Problem Theory: Methods for Data Fitting and Model Parameter 
Estimation," pp. 167-185 (1987). A Monte Carlo analysis begins with a 
computational model. In the model of present consideration, the accuracy of the input 
parameters is described in terms of probability distributions. In the Monte Carlo 



Express Mail No.: ET608077262US 



Date Filed: 01/09/02 



ATTORNEY DOCKET NO.: PM 98.086 



- 16- 

analysis, a single trial consists of randomly selecting one value from the probability 
distribution for each input parameter and calculating the resulting output. The 
complete simulation entails performing thousands of trials. For the Monte Carlo 
analysis specific to VLSA, the result of a given trial is retained only if that trial 
produces calculated logs that match the input log values within the specified closeness 
of fit. The results of all the retained trials are stored and analyzed to draw conclusions 
about the uncertainty in the solution of the original thin-bed problem. 
[0045] In the illustrative example, the accuracy of the input parameters for the 
interval-average VLSA analysis approach is shown in Figure 8. The parameters in the 
table of Figure 8 are from the model setup of Figure 7. If all the parameters of Figure 
7 were known exactly, then the base-case solution would provide the best and final 
answer for HPV. However, in practice, there is some uncertainty associated with each 
parameter. The entries tabulated in Figure 8 describe the uncertainty of each input 
parameter in the form of a bounded normal distribution. The numerical entries in the 
table include the mean values from Figure 7, the standard deviation (StdDev), 
minimum value (Min), maximum value (Max). The parameters Amin and Amax are 
calculated from Min = Mean - Amin x StdDev and Max = Mean + Amax x StdDev. 
The StdDev value determines how peaked or flat the distribution is between its Min 
and Max as illustrated in Figure 9. In general, a parameter that is known with high 
accuracy is characterized by a narrow, peaked distribution, while a parameter that is 
less constrained should be characterized by a broad, flat distribution. The values of 
Amin and Amax represent the number of standard deviations between the mean value 
and the min and max, respectively. As an example, Figure 9 shows three possible 
probability distributions for the total water saturation of high-quality sand parameter. 
[0046] For the example shown in Figure 8, the fractions of HiQ Sand, LoQ Sand, 
and Shale are all assigned uniform (flat) distributions in the range 0.00 - 0.70, 
representing a very high level of uncertainty in the values found in the base-case 
solution. The water saturation of both HiQ Sand and LoQ Sand are in uniform 
distributions of 0.05 - 0.30 and 0.25 - 0.45, respectively, again representing high 
uncertainty in these values. The aspect ratios are uniformly distributed in the range 
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0.0001 - 0.1000, indicating uncertainty in the degree of continuity of the thin sand and 
shale beds. Other parameters are restricted to relatively narrower ranges representing 
lower levels of uncertainty. 

[0047] Figure 10 shows a worksheet that gives a view of the values calculated for 
each Monte Carlo trial. For each trial, random samples are drawn from each of the 
distributions described in the table of Figure 8. With these parameters as inputs, the 
petrophysical forward model is re-calculated with results as shown in Figure 10. The 
joint probability index ("Joint PI") is calculated as a measure of the fit between the 
model-derived log averages and the input log averages. Depending on this measure of 
fit, the trial parameters are replicated to produce between zero and ten cases for later 
statistical analysis. The example shown was replicated to produce seven cases. The 
example shown was trial number 492 of 1000 "successful" trials - those trials whose 
measure of fit yielded one or more cases. 

[0048] After the Monte Carlo simulation, the dataset of retained cases is analyzed. 
Each of these retained cases represents a feasible solution to the thin-bed problem. 
The solution is feasible because it meets two conditions: (1) it is within the assumed 
distributions of the input parameters, and (2) it fits the measured log values to within 
the tolerance specified by the log standard deviations. The analysis generates and 
displays statistics on these feasible solutions. Figure 1 1 shows this statistical display 
for the combined HiQ and LoQ sand bed types. Counterclockwise from the upper 
left, statistics are presented for the hydrocarbon pore volume HPV(fit), the average 
sand porosity POR, the average sand water saturation Sw, and the total sand fraction, 
Vsd. In each plot, the shaded histograms represent the distribution of feasible 
solutions and the empty histograms represent the a priori distribution for the same 
parameter. 

[0049] Several observations about these statistics give an indication of the power 
of the VLSA method. First, the method recovered a good answer to the thin-bed 
problem that was posed: the central values (mean and median or P50) for each 
volumetric result are very close to the known interval-average values for the synthetic 
dataset, shown in Figure 6. Second, the method shows that the logs constrain the 
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solution for hydrocarbon pore volume to within a standard deviation of 5% of its 
value, despite the very high a priori uncertainties assumed for the bed-type fractions 
and water saturation. Likewise the logs constrain the solution for total sand fraction 
to within a standard deviation of 10% of its value. On the other hand, water saturation 
is not constrained at all relative to its a priori distribution, confirming the assertion 
that the conductivity log loses its sensitivity to water saturation in thinly bedded 
reservoirs. 

[0050] Although a particular, detailed embodiment of the present invention has 
been described herein, it should be understood that the invention is not restricted to 
the details of the preferred embodiment. It will be apparent to those skilled in the art 
that many modifications and variations to the embodiment described are possible 
without departing from the spirit and scope of the present invention. 
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